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1 Introduction 

In this work we study the time function r on the Schwarzschild black hole spacetime 
having the following properties: 

(i) The level sets of r result from evolution of a time- symmetric Cauchy slice of 
Schwarzschild by maximal surfaces under the additional requirement that proper 
time for asymptotic observers at infinity, which are at rest relative to the slicing, 
runs equally fast at both spatial ends. 

(ii) The time function r is zero on the time-symmetric slice and coincides with the 
proper time of the infinite observers. (This means that a, the lapse of the time 
function goes to one at both infinities along each slice.) 

Note that (i) is really a property only of the slicing defined by r rather than r itself. This 
time function which has first been considered in [jf], has two key properties: The first 
property is that r takes all real values or, in other words, the future singularity at r = 0 
does not prevent r from assuming arbitrarily large positive values (and similarly for the 
past). It is believed that this property holds on vacuum spacetimes more general than 
Schwarzschild. Here it is important to realize that such spacetimes are not “given” to 
us: rather they have to be generated by a Cauchy problem: One first constructs regular 
asymptotically flat initial data, satisfying the vacuum constraints, say maximal, and then 
tries to evolve these in time by analytical or numerical means. Doing this involves an a 
priori choice of gauge which in particular implies that the resultant globally hyperbolic 
spacetime comes already equipped with a specific time function. Suppose the initial data 
has a future-trapped surface. Then, by the Penrose singularity theorem 0, any Cauchy 
evolved spacetime is singular in the sense of having future-incomplete null geodesics. 
(Similar conclusions, but in both the future and the past direction, hold when the initial 
data has an outer-trapped surface [f|] or when the topology is nontrivial e.g. in the 
sense that there is more than one asymptotic end [5].) Many maximal initial data sets 
having one of these properties exist (for trapped surfaces, see |J). Now there is the 
conjecture, due to Moncrief and Eardley J7J, that if one evolves the initial data in a 
gauge where the whole slicing is maximal and r is proper time at infinity, the evolution 
should be extendable to arbitrarily large values of r, irrespective of whether singularities 
form or not. This global existence result, if true, would, in spirit at least, go a long way 
toward settling in the affirmative the Penrose Cosmic Censorshop hypotheses || in the 
case of asymptotically flat data. The spacetime evolving in the way described, in the 
Schwarzschild case, has the second property that it is in fact extendable: there are no 
maximal Cauchy slices of Schwarzschild reaching radii less than or equal to r = 3m/2. 
Thus maximal slices of Schwarzschild “avoid the singularity at r = 0”. It is this last 
property which numerical relativists expect to be true for evolutions of more general 
initial data and which is clearly desirable if numerical codes based on maximal slicings 
are used. 

Take any observer at rest relative to the slicing defined by r (“Eulerian observer”). 
Then / adr along the trajectory of that observer is her or his proper time. Since proper 
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time is finite as the slicing approaches the limiting maximal slice at r = 3m/ 2, we must 
have / adr < oo, and thus lim^oo a{r) = 0 (“collapse of the lapse” 0). Our main result 
is that, along the Eulcrian observers going through the bifurcation 2-sphere, 


a(r) 


4 A 


4 r 


3\/2 ° XP 3\/6 


exp 


3x76 


as r 


oo, 


m 


( 1 . 1 ) 


where the constant A is given by Equ. (3.41). The exponent in (1.1) has been estimated 
before [I], [TTJ by a mixture of numerical and model calculations. The estimate in [|Kj] 
of this exponent is 1.82 which agrees quite closely with our exact ~ 1.83. We hope 
that this result will be useful for the numerists as an accurate test for codes based on 
maximal slicings. An extension of the work here to the late time behaviour of a along 
the trajectories of arbitrary Eulerian observers will appear elsewhere [17 . 

Our plan is as follows: In § 2 we review some generalities on lapse functions and 
foliations. Then we give a precise definition of the time function under study. In § 3 
we perform the asymptotic analysis leading to Equ. (1.1). In Appendix A we essentially 
rederive the Schwarzschild metric in terms of spherically symmetric, maximal Cauchy 
data. In Appendix B we prove a calculus lemma which is basic for our analysis. 


2 Generalities 

Let (M, ds 2 ) be a globally hyperbolic spacetime and r : M —> R a time function, i.e. a 
function the level sets of which form a foliation T r of M by Cauchy surfaces = £. Then 
the function a : M —> R defined by 

a := (-(Vr) 2 )- 1 ' 2 (2.1) 

is called the lapse of T r . The reason for this name is that a measures “lapse of proper 
time” along trajectories normal to the leaves of T r as a function of r. To make this 
explicit, define the vector held by 

= -a 2 V M r => r^V^r = 1 (2.2) 

which is timclike and future (i.e. increasing r) pointing. The vector yields an orthogonal 
decomposition of M as M = R x £, as follows: Construct a diffeomorphism ip : R x E, 
i.e. <p : (A, y l ) G R x E i—► = <P\(y l ) 6 M, by 

<P\(y) = ^/P\ (y) = ^(ipx(y)) 

p 0 = t-\ 0). (2.3) 

Thus A, viewed as a function on M, coincides with r. We will, by abuse of notation, use 
the same letter r for A viewed in this way. Furthermore, since is normal to the leaves 
of JC T) the lines of constant y l are orthogonal trajectories, in other words 


, fir(y)Vr,i(v)S^(Vr(y)) = 0 . 


( 2 . 4 ) 









3 


Consequently, in the (r, ^-coordinates, the metric takes the form 

<p* T (ds 2 ) = (ffftg^dT 2 + <Pr, i <PTjgi*'dy l dy 3 

= 9 tt(t, y)dr 2 + ^(r, y)dy l dy\ (2.5) 

where y %] is the induced metric on the leaves and 

9rr(r,y) = -a 2 (ip T (y)). (2.6) 

Thus, along y l = constant, the proper time s is given by 

s = J a((f T t(y))dr'. (2.7) 

Note that, when t' is another time function giving the same foliation, i.e. t' = the 

lapse a changes according to a! = f ] a. Suppose now we are given another vector 

\drj 

field on M. This can be uniquely decomposed 

^ = Nn» + X 1 *, = 0 (2.8) 

where n 11 = — aV^r, is the future normal of jF r . To distinguish N from a, we call N the 

boost function of relative to T t - If X is non-zero on some leaf S TQ , it can be viewed 

as the restriction to E ro of the lapse of the time function t! obtained by = 0, 

t'lv = const. 

|Zjt 0 

We have the relation 

N = a fV fl r, (2.9) 

which is of course trivial in the present context, but will be extremely useful in our com¬ 
putation of the lapse a of a maximal foliation of the extended Schwarzschild spacetime, 
where £ ,J ' can be chosen as the ”static” Killing vector. 

We now recall some features of Schwarzschild which are used in our construction. In 
the exterior region r > 2m > 0 we have 

ds 2 = — ^1 — ‘—— S j dt 2 + ^1 — ‘—— S j dr 2 + r 2 dfl 2 , — oo < t < oo. (2.10) 

ds 2 can be smoothly extended across r = 2m to the Kruskal spacetime M on which r is 
a globally defined function r : M —> R + , which has saddle points at S, the bifurcation 
2-sphere of the horizon. The Killing vector field d/dt extends to a global Killing vector 
field on M which is spacclike in the interior, i.e. black and white hole, regions, null 
on the horizon and zero on S. Both the black hole region and the right exterior region 
can be written in the form (2.10) with the understanding that the functions ( 9 , tp) and r 
together with the retarded Eddington-Finkelstein coordinate 


u — t — r — 2m lg | r — 2m\ 


( 2 . 11 ) 
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covers both regions and the horizon at r = 2m. The function t goes to oo at the right 
component (where “right” refers to the original unextended spacetime) and goes to — oo 
at the left horizon. The set where t vanishes is the union of S, the original t = O-spacclikc 
hypersurface (extended in the obvious way to the left exterior region) and the timelike, 
totally geodesic cylinder T, which is ruled by timelike radial geodesics through S which 
are orthogonal to S, and which hit the singularity as r —» 0. Since r is constant along the 
trajectories of and r is, by (2.10), an “areal radius”, it follows that every spherically 
symmetric spacelike slice has a spherical minimal surface (a “throat”) exactly where it 
is tangential to £ /J ' (which of course can only happen in the interior and it necessarily has 
to happen there for slices leaving to the other (left) exterior region). 

Consider the function h(r, C ) given by 

k ^ C ^ = ~ l c (1 - 2m/x)(x 4 - 2mx 3 + C 2 ) 1 ^’ ( 2 ‘ 12 ) 

where the integral is to be understood in the Cauchy-principal-value sense for r > 2m 
and where 0 < C < 3-^ m 2 , r > rc and rc is the unique root of P(x ) = x 4 — 2mx 3 + C 2 
for this range of C in the interval 3m/2 < rc < 2m. For x > rc , we have P(x) > 0. Thus 
h(r , C) + r + 2m lg |r — 2m\ depends smoothly on (r, rc). We easily infer that 

t = h(r,C) (2.13) 

defines, for each fixed (7, a spacelike slice Sc which smoothly extends to the black hole, 
where it intersects T at r = re¬ 
in order to see that this surface extends smoothly and symmetrically through T, we 
use for r < 2m the parameter 


r = 


f r x 2 dx 
Ire [P(:r)]V2’ 


(2.14) 


which is the proper distance along the slice, as can either be seen from (2.10,12,13) or 
from Appendix A. Then, from (2.12,14) we have the system of ODE’s 


dh C 

dd r 2 — 2 mr 

d 2 r m 2 C 2 

d£2 r 2 r 5 ’ 


(2.15) 


with h( 0) = 0, r(0) = rc, ^(0) = 0, which is regular at i = 0. Thus the function 
r along the slice is symmetric with respect to t = 0 and smooth. This implies that 
dr/d£ = (1 — 2 m/r + C 2 /r^) 1 / 2 is antisymmetric. 

Next we observe that the level sets of o = t — h(r , C ), for fixed C in the allowed range, 
give rise to maximal surfaces on the Kruskal manifold, i.e. they satisfy 


W([-(Vct) 2 ]- 1 / 2 V„ct) = 0 . 


( 2 . 16 ) 
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The function a is not the time function of interest to us (in fact: cr being not differentiable 
at rc it does not define a global foliation). Rather this local foliation arises from moving 
a given maximal slice, say a = 0, along the flow of = ( d/dt ) M . The function N = 
[— (Vcr) 2 ] -1 / 2 is nothing but the boost function of d/dt relative to a = 0. In fact, in the 
explicit solution of (2.14) due to Reinhart |2J, he first, essentially by guessing, finds N to 
be 

V2\ 1/2 


N = 1 


2m C 2 
r r 4 , 


(2.17) 


and from this (2.12) can be inferred. For a more illustrative derivation from the initial- 
value point of view see Appendix A. Note that N as a function of l is antisymmetric 
relative to i = 0. 

We now claim that the surface t = h(r, C') lies everywhere in the future of t = h(r, C) 
when C > C ; and that t = h(r, C ) lies to the future of S' = S 0 . Amusingly we are unable 


to see this from the explicit integral (2.12). Instead, one first computes /j^rc from 


r% — 2 mr% + C 2 = 0, 


to yield 


dr, 


c 


2 C 


< 0 . 


(2.18) 


(2.19) 


dC 4r)f(l — 3m/2rc) 

Thus the claimed behaviour is true at least along the throat. Next observe that our slices 
are asymptotically flat at both spatial ends and that too(C') = lim^oo h(r, C) exists. 
Suppose that h(R, C) = h(R, C) for some R > r c to the right of T. Then, by the 
symmetry w.r.t. T this would have to happen also to the left of T. Thus we would have 
a lens-shaped region spanned by two maximal slices. But this, by an elegant argument 
due to Brill and Flaherty [|ll|], is impossible, except if the two slices are identical, which 
they are not in our case. This argument continues to be valid for R = oo. Thus t(r,C ) 
monotonically increases with C for fixed r, and so does too(C)- 

It follows that the equation t = h(r, C ) can be solved for C to yield a smooth time 
function defined on the r < 3m/2-subset of the part of Kruskal lying in the future of 
the Cauchy slice S. C gives the foliation we are interested in, but it is not yet the time 
function we want: rather this is obtained from eliminating C in terms of r in the relation 


T = too(C). 


( 2 . 20 ) 


Suppose we had started with the Cauchy slice t — 0 which, being time-symmetric, is 
in particular maximal and evolve it into a maximal slicing by a lapse function a going to 
1 at both spatial ends. This is possible in a unique way (see JT2[|). Then the resultant time 
function is spherically symmetric and symmetric w.r.t. T, and so it has to coincide with 
the one obtained above. In particular it follows that our r can be smoothly extended to 
negative values of r which would have been very non-obvious from the explicit formula 
( 2 . 12 ). 
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We next compute the lapse function a of r. Using (2.9) this involves computing 


(fV.r )- 1 


dC_ dh_ 
dr dC 


Note that (2.21) blows up at r = rc but in such a way that 


a = K'V,r)-‘iV 


( 2 . 21 ) 


( 2 . 22 ) 


has a smooth limit as r —> rc, as it has to be. Using formula (B.12), there results 


( dr \ 1 1 

1 

i 

r x(x — 3 m)dx 

VdcJ 2 

r — 3m/2 J 

r c (x — 3m/2) 2 [x A — 2mx 3 + C 2 } 1 / 2 


with 

dr r°° x{x — 3 m)dx 

d.C Jr c (x — 3m/2) 2 [x 4 — 2mx 3 + C 2 ] 1 / 2 
Note that N and a are linearly independent radial solutions of 

(A - KijK^f = 0, 


(2.23) 


(2.24) 


(2.25) 


where N goes to 1 at the right infinity and to —1 at the left one whereas a goes to 1 at 
both ends. 

We are interested in studying a along the trajectories of Eulerian observers. This 
requires choosing a coordinate p = p{r, r) the level surfaces of which are timelike cylinders 
orthogonal to our slicing. (One such timelike cylinder is already known, namely T given 
by r — rc-) Such a coordinate can be found without any calculation. Recall that maximal 
slicings preserve spatial volumes along Eulerian observers. Thus a suitable coordinate will 
be the “volume radius” on each slice, defined by 


P (r,r) = 3 


x 4 dx 


Jr c(T) [x 4 — 2 mx 3 + C 2 } 1 ! 2 ’ 
using that the spatial metric on each slice has the form (see Appendix A) 


2m C r 


-l 


gijdx^x 1 = 1-1-- dr 2 + r 2 dQ 2 . 




In the coordinates (t,p,6,ip), the Schwarzschild metric for r > 3m/2 reads 


ds 2 = —a 2 dr 2 + (—Y dp 2 + r 2 dQ 2 
\rj 


(2.26) 


(2.27) 


(2.28) 


where r = r(p, r) is given implicitly by (2.26) and a by (2.23). (To check (2.28) explicitly 


one should first observe that C 


= n 2 PP- 

P dC 
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Note that, as C approaches ^27/16 m 2 , rc approaches the value 3m/2, since 
P(x) =x 4 - 2mx 3 + C 2 =(x-^y (* 2 + rnx + ^ + O (^x - . (2.29) 


Equ. (2.29) also shows that rc approaches a double root of P(x) as C —> ^27/16 m 2 . 
Thus, as one lets r tend to infinity for fixed p, the function r approaches 3m/2. In that 
sense the slices approach the limiting maximal slice at r = 3m/2. We are interested in 
estimating a in that limit. For simplicity we will confine ourselves to p = 0, i.e. the throat 
T. 


3 The late time analysis 

It is convenient to replace the parameter C by 5 defined by 

5 — rc — ~rr, Tq — 2mr% + C 2 = 0. 


c 


(3.1) 


As C ranges between 0 and 3 -/p m 2 , 5 ranges monotonically from m/2 to 0. Using the 
rescaled quantities 


we find that 


where 


We have the 


C = 


c 


T 


T = 


nr 


m 


s = ~, 

m 


f{S) = -C 


ydy 


3 / 2+5 (y — 2)(p 4 — 2 y 3 + U 2 ) 1 / 2 ’ 


Q\ 3/ 2 /I \ V 2 


(3.2) 

(3.3) 

(3.4) 


Lemma: 


t(< 5) = —In<5 d—In |18(3\/2 — 4)| — 2 In 


3^-5 


9^-22 


+ 0(5) 


3^ 


lnh + A + 0(<5) as 5 —> 0. 


(3.5) 


Proof: First note that 

d 

dC 


C 


(p 4 — 2p 3 + C 2 ) 1 / 2 


y 3 (y- 2 ) 


(p 4 — 2p 3 + C* 2 ) 3 / 2 ’ 


(3.6) 



Thus, from the mean value theorem, 


y 


C 


^27/16 


| y - 2 ^ (r/ 4 - 2 y 3 + C 2 ) 1 / 2 ( y 4 _ 2j/ 3 + 27/16) 1 / 2 

where we have used 

27/16 — C < V27 P. 


< 


V27S 2 y 4 


(r/ 4 - 2 y 3 + C 2 ) 3 / 2 


(3.7) 


(3.8) 


(3.7) is valid for y ^ 2 but, by continuity, also for y — 2. We will find it convenient to 
sometimes express C in terms of 5, using (3.4). Writing 


Q(s) = s 2 (s 2 + 4s + ^ - <5 2 (V + 45 + ^ 
Equ. (3.3) can, after substituting s — y — 3/2, be written as 


. , 3\ 3 > 2 fl ,\ 1/2 r°° (s + 3/2)* 

T_l + 2/ U ) h (l/2^s)[Q(s)] 1 / 2 ' 


(3.9) 


(3.10) 


It is elementary to see that, for s > 5, 

0 < ^(s 2 - 5 2 ) < Q(s) < ( s 2 - 5 2 ) 

which, using a/1 + x < 1 + x/2 for x > 0, implies 


9 


+ 2s(4 + s) 


(3.11) 


[Q(s)] 1/2 [|(s 2 -5 2 )]V2 - [|(s 2 -5 2 )]i/2- 


2s 

< 9 


(4 + s) 


(3.12) 


The estimate (3.7) now takes the form 
I s + 3/2 f(5 + 3 / 2) 3//2 (1 / 2 - 5) 1 ' 2 


^27/16 


1 /2-s 


[ Q (*)] 1/2 


[s 2 (s 2 + 4s + 9/2)]V2 

v/27 <5 2 (s + 3/2) 4 


< 


V^7<J 2 (s + 3/2) 4 


[Q(s )] 3/2 


“ [|(s 2 -5 2 )] 3 / 2 


(3.13) 


The inequalities (3.12,13) are the basic estimates we will be using. We now split the 
integration domain in (3.10) as follows: 


5 < s < JS/2 , y<5/2 < s < 00 


(3.14) 


and write 


r = Ti + T-2 


( 3 . 15 ) 
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accordingly. We furthermore define (0 < 8 < 1/2) 

•s/Vi s + 3/2 


= v27/16 


ds 


r 2 ° = a/27/16 


Is (1/2 — s)[|(s 2 — 5 2 )] 1 / 2 
/■°° s + 3/2 

yi+ (1/2 - s)[s 2 (s 2 + 4s + 9/2)] V 2 


ds. 


(3.16) 

(3.17) 


Equ. (3.17) is in the principal-value sense at s = 1/2. These integrals can be explicitly 
computed using the formulas (see e.g. 


dx 


x 2 Vx 2 — 5 2 
dx 


\Jx 2 — 8 2 
x8 2 


x > 5 > 0, 


x 


2 \Jx 2 — S 2 


= In \x + x 2 — 5 2 \, x > 5 > 0, 


(3.18) 

(3.19) 


dx 


1 | — 2 Jc(ax 2 + bx + c) + 2c + bx\ 

_ In _I_ 


= —^ In 


Using 


x\J ax 2 + bx + c \fc 
s + 3/2 1 


2x 


c > 0 . 


(3.20) 


s(l/2 — s) s s — 1/2 


, there results after straightforward manipulations 


B° = 


To = 


In 8 H — In \J8/ 2 + o(l) as 8 —» 0 


(3.21) 


3VQ 


In \J8/2 + 


M ln2 + M ln 


18 

9 In 

3a/3 — 5 

4 + 3a/2 

A ill 

9y/6 — 22 


+ o(l) 


as 8 —> 0. 

Next we have to estimate the remainders. We have 

Ari = 


rV^ 

C(8) 

++ 6 

Is 

lQ(s)Y /2 

[|(s 2 — d 2 )] 1 / 2 


ds. 


(3.22) 


(3.23) 


Using C(8) = 27/16 + 0(8 2 ) and (3.11,12), this has 

r^/f/2 sc [ s 


Next 


A To = 


IAril A const 


s + 3/2 


= 0 ( 8 1/2 ). 


I Vs/ 2 1/2-s 


Is \/ s 2 — A 2 

C(8) a/27/ 16 


[Q(s)] 1 / 2 [s 2 (s 2 + 4s + 9/2)] 1 / 2 


ds. 


By (3.13), this has a bound of the form 


At 2 | < const 8 2 


(s + 3/2 ) A 


ly/JJl [s 2 (s 2 + 4s + 9/2)] 3 / 2 


ds = const 8 2 I. 


(3.24) 


(3.25) 


(3.26) 
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The integral / in (3.26) can be further split as / = / 2 + I ' 2 , where 

ri f«-4-3/9'i 4 /-i ds 


A = 


s + 3/2) 4 , r 1 

ds < const ' 


^ [s 2 (s 2 + 4s + 9/2)] 3 / 2 
df 


y/Zfi (S 2 - (5 2 ) 3 / 2 


l r 1 / 6 
= const-^T 


Now 


r* = 


8 2 J^rs (/ 2 - 1) 3 / 2 

(s + 3/2) 4 


h [s 2 (s 2 + 4s + 9/2)] 3 / 2 
Thus Ar 2 = 0(8). Putting all this together implies 


= 0 ( T ). 

o, 


r oo 

ds < const / — < oo. 

Ji s 2 


r(<5) = 


3^ 


In5 + A + o(l) as <5 —> 0, 


(3.27) 


(3.28) 


(3.29) 


which is not quite good enough. From Equ. (B.12) in the limit that r goes to infinity and 

2C %=-«( s+ !) 2 ' < 33 °) 

we see that 


df 8(8 + 3/2) A 2 r°° (s + 3/2)(s - 3/2) 

ds ~ (1/2 - sy / 2 h s 2 [g(s)]V 2 


ds 


= \/3 <5 


r 1 s 2 — 9/4 ds 
Is 1? [g(s)] 1 / 2 

J can in turn be split as J = J° + A J, where 

r 1 s 2 — 9/4 ds 


+ 0(S) = V3 8 J + 0(8). (3.31) 


r = 


\ 2.( c 2 


S' -(5 2 )] 1 / 2 




df 


h 


+ 0 7 


FvT^T VA 


0(ln(5) - y 9/8 -j • 1. 


Finally, 


A J = 


f l s 2 - 9/4 

1 

1 

Is s 2 

[[QM ] 1/2 

[|(s 2 -d 2 )] 1/2 _ 


ds. 


Thus, using (3.12) 


IAJI < const [ — 
Js s z 


ds 


const 


V s 2 — 8 2 8 


Putting (3.31,32,33) together, there results 


dr 3\/6 1 - 

7 TT =--— j + 0(1) as 8 —> 0. 

dS 4 8 w 


(3.32) 


(3.33) 


(3.34) 


(3.35) 
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Integrating (3.35), we obtain 


f(S) = In 5 + A'+ 0(5), 


(3.36) 


for some constant A'. Comparing with (3.29) we infer A = A 1 and the proof of the 
estimate (3.5) is complete. 


From A = A' and (3.36) it is elementary to infer that 

4 


5 = exp 


3V6 


(f - A) + O 


exp 


3v/6 


as r —> oo. 


Using (3.30), (3.35) can be written as 


dr _ 3 1 /1\ 

dC ~ 4^2 U, ‘ 


(3.37) 


(3.38) 


We want to evaluate the lapse a of the time function r = mr along the central throat 
r = rc■ This, using (2.23) and (2.24), is given by 


dr 


- i 


am = 


2 mS \ dC 


(3.39) 


Using (3.36) this finally leads to 

4 - 


a r = 


3\/2 

4 


5 + 0{5 2 ) 


wi exp {m) 

for r —> oo. 

We sum up our results in the 


exp 


4r 


3 y/6 


o 


m 


exp 


8 r 


3^ 


m 


(3.40) 


Theorem: For the chosen maximal foliation, with the time function r coinciding with 
proper time at infinity and being zero on the time symmetric leaf S, the lapse along the 
central geodesics orthogonal to the leaves behaves, as a function of r, according to (3.40) 
with A given by 


A = In 118(3\/2 — 4) | - 2 In 


3^-5 


9VQ-22 


= -0.2181. 


(3.41) 


It would be interesting to estimate the lapse for large r along arbitrary Eulerian observers 
rather than just the ones along T. In terms of the coordinate p introduced in § 2 we 
conjecture that 


= B{p) exp 


3^ 


+ 0 


m 


B 2 (p) exp 


8 r 

3^ m 


(3.42) 
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where B(p) behaves for large p as 


B(p) const • cosh 


4 /_p 
3^ 


(3.43) 


The form of B(p) in Eqn. (3.43) is motivated by the solution to the lapse equation (2.25) 
on the limiting slice at r = 3m/2, which is symmetrical with respect to the throat. 


Appendix A 


The following discussion is similar in spirit to |T5|. Let £ be the manifold R x S 2 with a 
Riemannian, spherically symmetric metric, which we write in the “radial” gauge, i.e. 


g = d£ 2 A r 2 (£)dQ 2 , r G (0, oo). 
The unit vector £ 1 = (d/dt} 1 is geodesic and satisfies (r' = drfdi) 

r' 

DiCj Qiji 


(A.l) 


(A.2) 


easily follows the lines of Besse [131, we find for the Riemann tensor 


where q t3 = g %3 —lii 3 and prime means derivative w.r.t. t. After a calculation, which most 

(A.3) 


Ri m lH k = - q u 


and 


Qi Qj Qk Qe Bi'j'k't' = -w(l - r r2 )q k[i q j]e . 


Identities (A.3,4) imply that 


„/2 


T \ — T 

Rij (2 £i£j A qij) A 3 Qij 

r " 1 _ r '2 
R = —4-b 2-—. 

rjn 


(A.4) 

(A.5) 

(A.6) 


The extrinsic curvature on £, in order to be spherically symmetric, has to be of the form 

Kij = v £i£j A w q i:j . (A.7) 

The condition K l3 g l] = 0 implies that v A 2 w = 0. Using (A.2) we have 


D % K %3 = ( v' A 3 — v 




(A.8) 
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Thus the maximal momentum constraint implies v = 2 C/r 3 for some constant C. Con¬ 
sequently, 


k- = « 

r 3 r 3 

C 2 

I\ ; j K '' = 6 —. 


(A.9) 

(A.10) 


Inserting (A. 10) and (A.6) into the Hamiltonian constraint, there results 


// -i 

r 1 — r 


/2 


— 4 — T 2 


Next we define m{r ) by 


„'2\ 


Now (A. 11) implies that dm/dr is zero. Thus 


6 7T- 

(A.11) 

C 2 

(A.12) 

2 r 3 

> V2 



r “( 1 _ T + u) ■ 


(A.13) 


Assuming m > 0 and 0 < \C\ < 3 m 2 , there are two initial-data sets consistent with 

(A.9) and (A. 13). One starts at r = 0, expands to an r max < 3m/2 and collapses back to 
r = 0. The other is an asymptotically flat complete metric on R x S 2 with mass m at 
both ends which is symmetric with respect to the throat at r = rc > 3m/2 with 


1 - 


2 m C 2 
- 1 —T 

r c 


0. 


(A.14) 


Here we restrict ourselves to asymptotically flat data. These constitute a 2-parameter 
family of solutions to the spherically symmetric, maximal vacuum constraints. Of course, 
we know from the Birkhoff theorem that members of this family with different C but 
the same m have all to lie in the same spacetime, namely the extended Schwarzschild 
spacetime. “Discovering” this fact in the present context amounts to finding the “height 
function” written down in Sect. 2. The trick is to try to find the remaining Killing vector 
and to seek the S^’s as graphs over the surfaces orthogonal to this Killing vector. If 
( gjj,Kij ) evolve to a spacetime having another Killing vector, there must be a function 
N , not identically zero, and a vector field X 1 so that 

2 NKij + 2D(,Xf ) = 0. (A. 15) 


Assuming X 1 to be again spherical, i.e. 

X\ )l f), 


d = d(r) 


(A.16) 
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and again using (A.2) and (A.9) we infer that 


2^ + 2r <£ 


NC 
4 —— + 2 r 


r 
, dfj, 


r° dr 

After combining (A.17) and (A.18), there results 


0 

0 . 


/z(r) 

N 


D 

^2 ’ 

D 

C 


r 


D = const 


(A.17) 
(A.18) 


(A.19) 
(A.20) 


where we have assumed C ^ 0. We assume without loss that D = C. The existence of 
(N,X t ) solving Equ.’s (A. 15) does not necessarily imply that the vacuum evolution of 
the initial data set has a static Killing vector. There also has to be satisfied 


C x K tj + DiDjN = N(Rij - 2 K u Kj). (A.21) 

It is straightforward to check that (A. 19,20) do satisfy (A.21). (In the case where C is 
zero, X 1 = 0 and Eq.(A.21) implies that N ~ r'.) 

We remark in passing that the function N, by virtue of (A.15) and (A.21), satisfies 


AN = NK ;j K ij . 


(A.22) 


(Of the two linearly independent spherical solutions of (A.22) N is that combination 
which vanishes on the throat.) 

It now follows that for r > rc the metrics 


ds 2 = -(N 2 - g ij X i X : ’)da 2 + 2g ij X j dx i da + g ij dx i dx j , (A.23) 


with Ni, X *, gij extended in a cr-independent way toRxE, are vacuum solutions evolving 
from the above initial data sets. They have = ( d/da) fl as a Killing vector. More 
explicitly, since 

2?77 

N 2 - XjX 1 = 1-, (A.24) 

r 

we have 

ds 2 = — ^1 — da 2 + 2 ^ dida + di 2 + r 2 d0 2 , (A.25) 

where r(£) is given implicitly by 


£(r) 


dx 


r c J1 — 2 m/x + C 2 /x A 


(A.26) 


(For C 


0 , £(r) can be written as £{r) = r^J 1 — 2 m/r + min 


l+^/l—2 m/r 
1— ^/l— 2m/r 
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Note that for C ^ 0 the above metrics extend smoothly across r = 2m. We now seek 
a function t with level surfaces orthogonal to d/da. Writing this function as 


t = F(r) + a, 


(A.27) 


we obtain from 


g^dx v = —(IV 2 — XiX^da + X^ = u(dF + da), 
for some function u, the equation 

- DF = ^k _ 

' JV 2 - XjXi 

which makes only sense off the horizon. Using (A.16,19) this leads to 

d£ _ C _1_ 

dr r 2 - 2mr ^1-2m/r + C 2 /r 4 ' 

Now consider the coordinate transformation 


(A.28) 


(A.29) 


(A.30) 


a — t — F. 


Then 

with 


ds 2 = -(N 2 - XjX j )dt 2 + ,j ij d.r i d.r J . 


9ij = 9ij + 2X(iF tj) — (N 2 — XgX^FjFj 


= gij + iN^^XeX^-'XiXj, 
where X, := gijXL Using (A.29) and (A.30), 


gijdx’ l dx J = 


1+1 


2m\- 1 C 2 




2 m\ 


-l 


dt + r 2 dn 2 = ( 1- dr z + r 2 dn 2 . 

r J 


(A.31) 
(A.32) 

(A.33) 
(A.34) 

(A.35) 


We have thus recovered the Schwarzschild metric. In particular this calculation shows 
that the parameter C in our initial-data sets is “pure gauge”: initial data with the same 
m lie in the same spacetime, namely as level sets of the function a. They can also be 
written as 

t = F c (r) (A.36) 

and its translates under = ( d/dt ) /i , where 


F c (r ) 


-C 



dx 

(1 — m/x)( 1 — 2 m/x + C 2 /x^) 1 / 2 


(A.37) 


for some r 0 . Taking r 0 = Tq we have, with h(r, C) = Fair), recovered (2.12). 
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It is shown in Sect. 2 that 

t = h(r , C) (A.38) 

implicitly defines a smooth time function on the r > 3m/2-subset of the future half of 
Kruskal. The boost function N obtained in this Appendix satisfies the same equation, on 
each leaf Sc as the lapse function of C, namely (A.22). The reason for this is that, for 
fixed Sc, defines another local foliation, which is again maximal since is a Killing 
vector. 


Appendix B 

Consider 


F(x,E) = 


W(y) 


;dy , 


(B.l) 


I xe =v-i(e) [E - V(y)]V 2 
where V is a smooth function V : [xo, oo) —> R with 

0 < V(x 0 ), V'(x) < 0 for x > x 0 , V(x) =0 (B.2) 

and 

0 < E < V(xo). (B.3) 

The function W is smooth except perhaps at x — x, where it may have a simple pole. 
(Thus the pole of \[E W(y)/[E — V(y)}^ 2 is independent of E.) In the latter case (B.l) 
is to be understood in the principal-value sense and the following operations valid for 
x 7 ^ x. Next define (we follow [l 6 j in spirit) 


J(x, E) = He- V(y)] 1 / 2 V(y)W(y)dy. 
Jx E 


Note that V ■ W is smooth. Equ. (B.4) can be rewritten as follows 


dy 

Differentiating (B. 6 ) w.r.t. E twice, we obtain 

1 


V'{y) 


dy , 


\E - V(y)fV 


d 

dy 


V(y)W(y) 

V'(y) 


d 2 1 

J(x,E) = ---. 


dE 2 


V(x)W(x) 1 


f x 1 d 

\V(y)W(y) 1 

L [E — V (y )] 1 / 2 dy 

L V{y) \ 


(B.4) 

(B.5) 
dy. (B. 6 ) 

dy. 


On the other hand, differentiating (B.4) once w.r.t. E it follows that 

V(y) 


1E J ^ = \ 


*b [B - V( y)]‘/2 

1 r x V(y) - E + E 

2 J XE [E-V(y)]V 2 


W(y)dy 
W (y)dy 


1 

2 


'X E 


[B - V(y)]V 2 W(y)dy+-EF(x,E). 


(B. 8 ) 
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Differentiating (B.8) once more w.r.t. E and comparing with (B.7) we finally find 


l -F(x,E)+ 1 -E-^F(x,E) = 

1 1 V{x)W{x) 

~2\E -V(x)] 1 / 2 V'{x) 


1 r x 1 d 

2 J XE [E — V (y)] 1 / 2 ~dy 


V{y)W{y) 

V'(y) 


dy. (B.9) 


In our case we will have that V'{xq) = 0 and we study the blow-up of F 0O (£') = 
lim^oo F(x, E) as E tends to E 0 = V(x 0 ). As for a mechanical analogue, we could 
think of a particle on a half-line in a repulsive potential V (x) and imagine F(x, E) to 
be the time it takes a particle of energy E to travel from Xo to x. (If it were not for the 
presence of W(y) in (B.l), this interpretation would be literally true.) The force on the 
particle grows so fast for large x , that the particle reaches infinity in finite time F oc (E). 
There is an unstable equilibrium point at x — xq. We ask for the way in which F 00 (E) 
blows up as E approaches V(xo). If the energy E is further increased,the orbits reach 
x = 0: this corresponds to maximal slices hitting the singularity. 

To make contact with our function h(r,C), set 


V{x) = -x 4 + 2 mi 3 , E = C 2 , W(x) = 

h(r ., C) = CF{r , C 2 ), 

3 m _ 

x 0 = —, x = 2m. 


1 — 2 m jx' 


(B.10) 


Thus 


± h ^ c ) = 2E± F ir, E ) 


+ F(r, C 2 ) 


(B.ll) 


E=C 2 


which, combined with (B.9), gives 


d 

dC 


h(r, C ) 


1 1 r r x(x — 3 m)dx 

2 (r - 3m/2)^/l - 2 m/r + C 2 /r 4 2 Jr c (x - 3m/2) 2 (x 4 - 2mx 3 + C 2 ) 1 / 2 ' 

(B.12) 
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